Statistics of a hydrophobic chain near a hydrophobic boundary 
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We study the behaviour of a hydrophobic chain near a hydrophobic boundary in two dimensions, 
using the decorated lattice model of Berkema and Widom [G.T. Barkema and B. Widom, J. Chem. 
Phys. 113, 2349 (2000)] to obtain effective, temperature dependent intrachain and chain-boundary 
interactions. We use these interactions to construct two model hamiltonians which can be solved 
exactly. Our results compare favorably with preliminary Monte Carlo computations, using the 
same effective interactions. At relatively low temperatures and at high temperatures, we find that 
the chain is randomly configured in the ambient water, and detached from the wall, whereas at 
intermediate temperatures it adsorbs onto the wall in a stretched or partially folded state, again 
depending upon the temperature, and the energy of solvation. 

PACS numbers: 5.50.+q, 64.75.+g, 87.15.Aa 



I. INTRODUCTION 

Non-polar molecules placed in water behave as if there 
were attractive interactions between them, at least within 
a certain temperature interval. This is the co called 
hydrophobic interaction and is entropy driven. |Q, |^] 
Hydrophobic interactions play an important role in bi- 
ological processes, most prominently in protein fold- 
ing 1 1, 1^, 1^ . Although many short proteins are able fold 
on their own within very short times, others require the 
help of chaperons to be able to do so. To be able to elu- 
cidate the role of chaperons, we believe that it would be 
useful to understand the behavior of a hydrophobic chain 
in the presence of a hydrophobic boundary. 

In trying to understand the behaviour of a hydrophobic 
chain in water, one must take into account both the 
hydrophobic interactions mediated by the orientational 
entropy of the water molecules, and the configura- 
tional entropy of the chain, while respecting its con- 
nectivity. To compute the temperature dependent ef- 
fective interaction potentials between the hydrophobic 
wall and hydrophobic polymer chain and also between 
the monomers of the hydrophobic chain, we borrowed the 
decorated lattice model recently proposed by Widom and 
coworkers [|[ 0, || . We will calculate the effective free en- 
ergy cost of bringing a hydrophobic solute molecule from 
the bulk to a distance r from the hydrophobic wall within 
the decorated lattice model in one dimension. This will 
provide an estimate for the interaction potential between 
the hydrophobic residues of the polymer chain with the 
hydrophobic wall in two dimensions. We will use the 
Mean Field Approximation(MFA) to find the effective 
interaction between neighboring hydrophobic residues in 
arbitrary dimension. 

To be able to compute the partition function of the 
hydrophobic chain in the presence of a boundary in two 
dimensions, we introduce two exactly solvable models, 
using these effective interactions. The first is a solid 
on solid (SOS) type of simplification of the allowed con- 
figurations of the chain. This allows us to reduce the 



problem to a one-dimensional lattice model, and use the 
transfer matrix formalism. The second simplified model 
again consists of performing exact sums over subsets of 
allowed configurations to obtain an approximation to the 
true partition function. Thirdly, we will perform Monte 
Carlo simulations. 

To investigate the thermodynamic properties of the 
chain we have focussed on the mean distance (rem) from 
the wall, and the average length (L) of the projection of 
the chain on the boundary. 

The paper is organized as follows. In section 2, we 
outline the lattice model within which the effective hy- 
drophobic interactions are computed. In section 3, we 
will discuss the various simplified models within which 
we have performed the sums over the chain configura- 
tions. In section 4 we discuss our results. 



II. MODELLING HYDROPHOBIC 
INTERACTIONS 

A decorated lattice model that mimics the solvent me- 
diated hydrophobic interaction was suggested by Widom 
and his collaborators 0, ^. In this model, g-state 
Potts spins, {si}, are situated at lattice sites. These rep- 
resent the polar solvent molecules. They can have any 
of the q different polarization directions. Hydrophobic 
molecules(HM) can only be accommodated at interstitial 
sites, more precisely on the bonds connecting neighboring 
pairs. Lattice-gas variables, {(Jij}, (Jij = (0,1), located 
on the bonds (ij), indicate whether an interstitial site is 
empty or occupied by a HM. 

Interaction between water molecules and HM is always 
attractive, because of the dipole-induced dipole interac- 
tion. On the other hand water molecules can form short 
lived tetrahedral structures [0 |ll| stabilized by hydro- 
gen bonds i.e., a type of short ranged order. Be- 
cause these structures have an open cage like space be- 
tween them, HM can be accommodated there with- 
out breaking any hydrogen bonds. Thus, this "ordered" 
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configuration is the minimum energy configuration of wa- 
ter molecules in the presence of HM. If there are no HM 
between the ordered water molecules, there still is an at- 
tractive interaction due to the hydrogen bonds and the 
dipole-dipole interactions, but the absolute value of the 
interaction energy is smaller, by precisely the amount 
contributed by the induced dipole interactions. At higher 
temperatures, water molecules will tend to be oriented 
randomly. This state, with no HM intermixed with the 
water molecules, is chosen as the reference, i.e., the zero 
level of the energy. When water molecules are randomly 
oriented, they can still have hydrogen bonds between 
them, though fewer in comparison to the ordered state. 
However, unlike the ordered state, there will be less open 
space between them. To be able to accommodate a HM 
in a disordered region of water molecules, further hydro- 
gen bonds have to be broken. Thus, the insertion of HM 
within this disordered phase of water molecules is ener- 
getically unfavorable. The Hamiltonian for the water- 
hydrophobic solute system can be written as M , 



<ij> 

aij w (l - Ss^^sjSs^^i)] ■ 
The interaction energies are ordered thus, 

w < u < < V , 



(1) 
(2) 

where v may be thought of as the solvation energy of the 
HM in the disordered state of the water molecules. In 
this model, the unique ordered state of the tetrahedrally 
bonded pentamers ]l3| , which is able to accommodate the 
HM without breaking any hydrogen bonds, is identified 
with the configuration where all the Si are in the state 1. 

We immediately realize that Eq.(^ may be rewritten 
in terms of two-state variables ti, defined by 

Ss,s = U. (3) 

where U — {1,0}. In the partition function the mul- 
tiplicity of the Si ^ 1 states can be taken care of by 
inserting a factor of {q — 1) for each Potts spin not in 
the ordered state, or a term — ti) \n(q — 1) into the 

Hamiltonian, to get, in one dimension, 

N 

H = {tiU+i [a.i{w ~u-v)+v]+aiV 

i 

-p-\l~U)\n{q-l)} (4) 
N 

= ''^^ Hi[ti,ti+l,<7i] . 
i 

(5) 



A. Effect of the boundary in 1-d 

Let us first consider a one dimensional system, in or- 
der to be able to estimate in closed form the effective 
interaction of a HM with the hydrophobic boundary. 



For N being the length of the one-dimensional lattice of 
water molecules, the free energy cost F{N, T, r) of adding 
only one HM at an interstitial a distance r from the wall 
at temperature T is given by 

where f3~^ = ksT as usual, Zq{N,T) is the partition 
function of the one dimensional system with cr^ = for 
all i, that is, no HM molecules, and Z{N, T, r) is the par- 
tition function computed in the presence of one HM a 
distance r from the wall. The effective interaction be- 
tween the wall and a single HM is thus given by the free 
energy cost of bringing HM from bulk to distance r from 
the wall, 

F^'^ (1, r) - FiN, T, r) - F{N, T, n) , (7) 

where means a displacement from the wall beyond 
which the effect of the wall is no longer perceptible, 
namely a bulk site. In the thermodynamic limit 



hm F'^>il,r) 

N^OQ 



(8) 



To compute the partition functions in (^, we used the 
transfer matrix method. From Eq.(^, the transfer ma- 
trices in one dimension are obtained as 



(9) 



Thus, the transfer matrix is conditional on the presence 
(or absence) of an interstitial HM at each bond connect- 
ing two water molecules, and we find, 



r(o) 
r(i) 



-/3 m 



(9-1)^ 
(g-l)i {q-1) 



-13 V 



(9-1) 
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-I3v 



(9-1) 



(10) 
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for the two possible resulting transfer matrices. ^From 
Eq. (0), we get (with one HM inserted at a distance r 
from the wall), 

-PFj^\l,r) = ln^(l|r'^-i(0)r(l)T(0)^-nfc) 

k 

-ln^(l|r^-i(0)r(l)|m). (12) 

m 

Notice that the left-most vector is fixed to be unity, sig- 
nalling the presence of the hydrophobic wall. In the ther- 
modynamic limit N oo, this reduces to, 



PF^^'>{l,r) 



ln^AiTij{l)ajiaik 

ijk 



In Y]aiiaijTji(l) , (13) 
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where we have defined 



= aiiau + {X2/Xiy ^021024 , (14) 



with 



A1.2 = -e-'^"(l + (g-l)e''" 

± [l-{q~l){2e~f'^~{q + 3)e^^nV) (15) 

being the eigenvalues of T(0), and a^i the elements of its 
fcth eigenvector. 

The free energy cost of bringing a pair of HMs from 
the bulk to a distance r from the hydrophobic wall could 
also be found, similarly, from taking the thermodynamic 
limit in 



F(^)(2,r)= hm F^^\2,r) , 



N 



(16) 



where, 



k 

-ln^(l|r^-2T2(l)|m), (17) 

m 

and T(0) is to be understood where the argument has 
been omitted. 

We have plotted F^-^^{l,r) in Fig. (2) as a function 
of the inverse temperature. We find that, at intermedi- 
ate temperatures F{l,r) becomes attractive. Although 
F^-^'>(2,r) is always positive, it becomes smaller at inter- 
mediate temperatures. Hydrophobic interactions become 
effective in a temperature interval which can be seen to 
range approximately between \u\/(2kB) < T < 2\u\/kB- 

We will use F^-^^{l,r), which we have calculated ex- 
actly in one dimension, to give us an estimate of the 
interaction between the HM and the hydrophobic bound- 
ary in two dimensions. On the other hand, to take into 
account the self-interactions of the chain in two dimen- 
sions we also need an effective temperature dependent 
pair potential. This is the subject of the next section. 



B. EflFective hydrophobic pair interaction in the 
Mean Field Approximation 

We made use of the Mean Field Approximation (MFA) 
to find the solvent mediated interaction between the hy- 
drophobic molecules on a cubic latice in arbitrary di- 
mension. In the MFA, the Hamiltonian (|l]) of a wa- 
ter molecule interacting with neighboring hydrophobes 
in two dimensions, in the presence of a mean field (t) due 
to the ordering of the ambient water molecules, can be 
written as, 

2d 

F[mf = ^ {t{t) [fTj {w — u - v) +u]+ a.j v} 

-r V(i - M-? - 1) . (18) 



where d = 2 and /j, has been inserted for later conve- 
nience in computing expectation values, and will be set 
to unity otherwise. Performing the t sum in the parti- 
tion function leads to couplings between the aj situated 
on the bonds emanating from the lattice site on which 
t is located. We define the effective two-body interac- 
tion strength M between the nearest neighbor (nn) and 
next nearest neighbor (nnn) pairs (which are indistin- 
guishable from each other in this approximation), as well 
as the four-body (or plaquette) interactions between the 
HM, via. 



(19) 



where [ij) denotes nn and nnn pairs. To the first ap- 
proximation 1^ we will neglect the plaquette coupling 
K. Keeping only the two-body interactions, we find. 



From 



^-2l3[{t){w~v+u)+v] _^ (,j-l)e-2/3" 
|-g-/3[(t)(t«-t.+3 «)+!)] ^_ _ X)g-/3i,j2 

(e-4^"<*> + g-l) . (20) 



= ^ln^U=i (21) 

one gets a self-consistent expression for (i), which we 
solved numerically for each given temperature T . Substi- 
tuting in Eq.(|2C|) one finally obtains the effective solute- 
solute inter actio n energy in two dimensions, which we 
plot in Fig.( |ll A| ) for different choices of q. The interac- 
tion between HMs is attractive for any finite tempera- 
ture. 



III. STATISTICAL MECHANICS OF POLYMER 
CHAIN 

We are interested in the behavior of a hydrophobic 
polymer chain in the presence of a hydrophobic wall. We 
used three different approaches. 



A. Modular chain or SOS model 

We define a set of elementary modules, from which a 
large number of chain conformations can be built, such 
that only nearest neighbor modules come within the in- 
teraction range of each other. The subset of configura- 
tions that can be generated by random combinations of 
the modules that are shown in Figure ( ^) can clearly 
be seen as graphs (taking the boundary as the axis) 
without overhangs, as in a restricted solid-on-solid (SOS) 
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model 11 in (1+1) dimensions, where successive steps 
are constrained to differ by at most one unit of height. 
Making use of the Hnearity of the chain and the restric- 
tion to nearest neighbor interactions between the mod- 
ules, we used the transfer matrix along the chain to solve 
the partition function exactly for our model Hamiltonian. 

We labeled the modules in Figure ( ^) as 1, 2, 3 from 
left to right. A chain configuration is uniquely specified 
by associating a variable, Ui = {1, 2, 3}, i = 1, . . . , Nm, 
with each module along the chain, and by specifying the 
distance of the first module from the wall. Note that the 
number of residues along the chain is given by 2Nm in 
this case. The interaction energy of each residue with 
hydrophobic wall is given in terms of F^-^'>{l,r). More- 
over, we took M{P), defined in Eg . (|20|) , to be the inter- 
action energy between nearest and next nearest neighbor 
residues, (see Fig.(^b)). Note that the nearest neigh- 
bor interaction (wavy line) connects residues belonging 
to modules twice removed from each other. Yet, since 
this occurs only in the + 1) = (2,3) or (3,2) combi- 
nation, independently of the identity of the i — 1st mod- 
ule, it can still be accomodated within a nearest-neighbor 
Hamiltonian. 

We model the effective Hamiltonian of a polymer with 
Nm modules as, 



i 

+ ^l(?'i-l,n) + /i2(n,'"j)} 



(22) 



The vectors \uii correspond to the three states of the 
variable u^, i.e., (1,0,0), (0,1,0) etc., so that the coeffi- 
cient of the pair interaction M{P) is conveniently written 
in terms of 



r = 




(23) 



The second term is the free energy cost of adding HMs to 
the solvent matrix, hi{ri) = — ^(-^^(1, r,j_i) - F(^'(l,ri). 
The distance of the second residue on the zth module 
from the wall, r^, is found from ^ r + pi, where r is 
the distance of the first module from the wall, and 



J2 ('^".^2 - <5„,-,3) 



(24) 



Note that the displacement of the first residue on the 
ith module is the same as that of the second residue on 
the i — 1st module, and therefore the expression for hi 
follows. The final term is 

h2{u^)^-F'i2,r, + l)Su^a-F\2,n-l)Su„3 , (25) 

and adds on the extra cost of placing a pair of HM in 
the same "row" perpendicular to the wall. In practice, 
we found that the addition of this term overestimates the 



interaction with the wall. While its inclusion or omission 
had a minimal effect on the average center of mass dis- 
placement from the wall, its inclusion led to unrealistic 
results for the projected chain length. Therefore in the 
calculations reported below, it has been set to zero. 
The partition fmiction of the polymer is. 



r {„} 



(26) 



Explicitly, 



Z = ^ {ri,ui\U\r2,U2){r2,U2\U\r3,U3) . . . 

• ■ • {rN„,^i,UN„,^i\U\rN^,UN^^) . (27) 

Here, \ri,Ui) are Ai x 3 dimensional vectors, with Ai 
being the size of the system in the direction orthoganal 
to the wall. The transfer matrix U is given by a direct 
product 



with 



(1) 



k£ 



U ■ 



3 

E 

C=i 



(28) 



(29) 



W, 



(2) _ 



5,,2[^^^''' [5,^1+ Ska) + e^^''5,,3] (30) 



W(3) = (2^3) 



(31) 



where k,£ — 1,2,3, (2 ^ 3) indicates an interchange of 
the indices 2 and 3 in the previous equation, and 

+ '5^3(7^^) , (32) 

where 7, 77 = 1, ... and C = 1, 2, 3. We note that only 
the diagonal, upper diagonal and lower diagonal elements 
of the matrices R^'^^ are different from zero. 
We calculated the center of mass distance, 



(33) 



of the hydrophobic polymer from the hydrophobic bound- 
ary as a function of temperature, and this is shown 
in Fig.(^). At intermediate temperatures the polymer 
chains are attracted to the wall so strongly that (rem) — 
1.5. The chains are predominantly in a zig-zag configu- 
ration confined very close to the wall, with half of them 
actually adsorbed on the wall, and the maximum number 
of nn and nnn interactions. 

As /3 ^ (high temperatures) the intrachain interac- 
tions M also go to zero, the entropy of the chain becomes 
the determining factor, and the chain floats free. At low 
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temperatures, as the entropy term in the free energy be- 
comes neghgible, the equihbrium state is determined by 
energetic considerations, and the polymers desorb and 
take on random configurations, constraining a large num- 
ber of water molecules in their neighborhood. 

The average end to end distance of the polymer chain, 
projected on to boundary, is given by 



(L) = N, 




(34) 



The temperature dependence is reported in Fig. (|6|). In 
the limit (3^0, clearly (L) = N^il + 1/3), which is 
what one sees in Fig.(||), with iV„ = 30. It is interest- 
ing to note the non-monotonic behaviour of (L) within 
the region of interest, namely the temperature interval 
for which the center of mass lies very close to the wall. 
This non-monotonicity arises from the competition be- 
tween the entropy mediated effective self-interaction of 
the chain (leading to smaller L) and the interaction with 
the wall (completely shielding one side from the water by 
stretching out to adsorb on to the wall) . This behaviour 
is also observed in the models we have considered in the 
subsequent sections. 

Although the SOS model is exactly solvable, it is un- 
able to take into account configurations of the chain 
which fold on themselves, and we therefore have also con- 
sidered a model where such conformations are allowed. 



B. The ri-fold model 

In this section we take a different tack, and perform 
exact summations over a richer conformation set for the 
same interactions, although this still includes only a frac- 
tion of all possible polymer configurations. These con- 
figurations are shown in Fig.(0). If the length of the 
polymer is Ni then the energy of a chain with an integer 
number of folds Ni /n, is given by 



Ni 



VF(^)(l,r + i) - Miynil-5ni-5nN,) (35) 



4=1 



where r is the distance from the wall (see Fig.(|^)) and f„ 
is the total number of nearest neighbor and next nearest 
neighbor pairs in this configuration, iy„ = 3(n — l)(7V;/n— 
1). We calculated (L) — {Ni/n), which is the mean 
value of the vertical distance between the first and last 
monomer, from 



(36) 



where the prime indicates that the summation is only 
over the subset of configurations described above, with 
Ni/n integer. We also calculated the center of mass dis- 
placement from the wall, (rem)- These quantities are 
shown in Figs.(||,||). 



In the low temperature limit, increasing the number 
of HM-water nn pairs lowers the energy and this is fa- 
vorable since the entropy term in the free energy is sup- 
pressed, and the chain takes on relatively open, random 
configurations. At intermediate temperatures where the 
hydrophobic interactions are the most effective, the chain 
prefers to neighbor the hydrophobic wall at as many nn 
sites as possible, and therefore is adsorbed on the wall in 
the unfolded state. As the temperature is raised some- 
what more, effective self interactions of the chain become 
more important, and the chain is in a more folded state, 
although still adhering close to the wall. At high temper- 
atures it is advantageous to minimize the number of near- 
est neighbor sites at which the chain is in contact with 
water molecules, since the entropy of the water molecules 
is rather large, especially for large q. On the other hand, 
the entropy of the chain also favors open configuration, 
which wins out in the high temperature limit. It should 
be noted that in Fig. (9), with Ni = 50, (L) is close 

to Nf^'^ — 18.8, at both extremes, with the power being 
that of the Self Avoiding Walk in two dimensions. 



Monte Carlo simulations 



For comparison, we display in Figs.(|0 11) preliminary 
Monte Carlo simulation results, for 3x10^^ random Self 
Avoiding Walk configurations of length N — 20. The 
initial points of the walks have been randomly chosen 
within a 100 x 100 square lattice. If a random walk 
passes through any lattice point which it has already vis- 
ited, the configuration is discarded, and a new one gen- 
erated. Each successfully generated configuration was 
decorated with the interaction potentials found in sec- 
tion (H), namely, F(^)(l,r) and M(/3) (Eqs.( 0,|o|)), to 
finally compute the expectation values for the center of 
mass displacement from the wall and the longitudinal 
component of the end to end distance, in the canonical 
ensemble. We will be reporting on more extensive Monte 
Carlo simulations in a separate publication. 



IV. DISCUSSION 

We have presented two exactly solvable versions of a 
model for the statistics of a hydrophobic polymer chain 
in water, in the presence of a hydrophobic boundary. The 
model of Widom and co-workers |^, ^ for hydrophobic 
interactions has been our point of departure for comput- 
ing approximate effective intrachain and chain-boundary 
potentials. Although the behaviour of chains (or mem- 
branes) in the vicinity of spatial boundaries have been 
considered before ||l^, |l6|, , these studies have concen- 
trated on temperature independent interactions. 

With the inclusion, to various degrees of accuracy, of 
the entropy of the chain, we are able to take into ac- 
count the competition between the entropy of the water 
molecules which can be constrained by the presence of 
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hydrophobic molecules in their neighborhood, and the 
entropy of the chain. We find that although at low and 
high temperatures, the chain prefers to be in a random 
configuration, detached from the wall, there is an inter- 
mediate temperature range where it is adsorbed on to 
the wall, at least for the relative values of the hydrogen 
bond, dipole-induced dipole and solvation energies which 
we have assumed. For relatively smaller values of the sol- 
vation energy, v, there is a sub-interval of temperatures, 
where the chain adheres to the wall in a relatively folded 
state. It is gratifying to find that this qualitative feature 
found in the exactly soluble SOS model and the n-fold 
model is reproduced in the Monte Carlo calculation, al- 
though the adsorbed "phase" occurs at somewhat higher 
temperatures. 



The present study provides one way of modelling the 
interaction of a protein with a hydrophobic surface. Our 
results are interesting from the point of view of pro- 
tein dynamics, especially for the folding of relatively long 
polymer chains, which may need the assistance of chap- 
erons to fold correctly. The effect of hydrophobic 
boundaries on the folding of amino-acid chains is also 
interesting from an evolutionary point of view, as has 
been suggested by Tiizel and Erzan |T^. Further work 
is in progress, to include polar as well as hydrophobic 
elements, for a more realistic protein chain. 
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FIG. 1: Decorated lattice model. Lattice sites are occu- 
pied by water molecules (shown as filled circles) ; hydrophobic 
molecules (open cirlces) can only be accommodated at inter- 
stitial sites. 
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FIG. 2: The effective interaction potential of a residue with 
the hydrophobic wall for different values of q, at r = 1, at dif- 
ferent inverse temperatures. Here q is the number of different 
orientations which can be assumed by the water molecules. 
The interaction coefficients of the lattice model were chosen 
to be i(j = —1.5, u = —1, V = 1. 



FIG. 3: Effective, temperature dependent nn and nnn inter- 
action energies between hydrophobic residues in water, in the 
MF approximation to the decorated lattice model [^. Differ- 
ent q values are shown. The effective interaction is stronger 
for larger q. The coupling constants for the decorated lattice 
model are the same as in Fig. 2. 



(1) 



(2) 



(3) 



r. 



FIG. 4: a) (Top panel) Elementary modules used to generate 
SOS like chain configurations which only allow nearest neigh- 
bor interactions between the modules, via nn or nnn inter- 
actions between the hydrophobic residues, b) (Lower panel) 
Nearest and next nearest neighbor interactions M{/3) between 
HMs on the chain are indicated as wavy and dashed lines, re- 
spectively. 
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FIG. 5: The average center of meiss displacement from the 
boundary, of the hydrophobic chain with 60 residues in the 
SOS approximation, for different values of solvation energy, 
V, and different values of q. For computational purposes, the 
width of the channel was chosen to be 12 lattice spacings. 
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FIG. 6; Mean length of the hydrophobic chain with 60 
residues, projected on to the boundary, in the SOS approxi- 
mation. Different values of solvation energy, v, and different 
values of q are shown for comparison. 
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FIG. 7: Polymer configurations included in the exact enumer- 
ation of the n-fold model. 



FIG. 8: The average center of mass displacement from the 
wall, of a polymer with 50 residues, for different values of the 
solvation energy, v, and of q, in the n-fold model. The width 
of the channel was taken to be 50 lattice units. 
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FIG. 9: The average longitudinal component of the length of 
the hydrophobic polymer with 50 residues, for different values 
of solvation energy, v, and of g, in the n-fold model. 



Plul 



FIG. 10: Preliminary Monte Carlo results for average center 
of mass displacement of the hydrophobic chain from the hy- 
drophobic boundary, for different values of solvation energy, 
V. The chain length is 20 residues, and the channel size is 100 
lattice spacings. 
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FIG. 11: Same as in Fig. 10, for the projection of the end to 

end distance of the hydrophobic chain on the boundary, for 
different values of the solvation energy, v. 



